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Spatially resolved study of the physical properties of the ionized gas 
in NGC 595 
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ABSTRACT 

We present Integral Field Spectroscopy (IFS) of NGC 595, one of the most luminous H II 
regions in M33. This type of observations allows us to study the variation of the principal 
emission-line ratios across the surface of the nebula. At each position of the field of view, we 
fit the main emission-line features of the spectrum within the spectral range 3650-6990 A, 
and create maps of the principal emission-line ratios for the total surface of the region. The 
extinction map derived from the Balmer decrement and the absorbed Ha luminosity show 
good spatial correlation with the 24 /im emission from Spitzer. We also show here the capa- 
bility of the IFS to study the existence of Wolf-Rayet (WR) stars, identifying the previously 
catalogued WR stars and detecting a new candidate towards the north of the region. The ion- 
ization structure of the region nicely follows the Ha shell morphology and is clearly related 
to the location of the central ionizing stars. The electron density distribution does not show 
strong variations within the H II region nor any trend with the Ha emission distribution. We 
study the behaviour within the H II region of several classical emission-line ratios proposed as 
metallicity calibrators: while [Nll]/Ha and [Nll]/[Olll] show important variations, the R23 
index is substantially constant across the surface of the nebula, despite the strong variation of 
the ionization parameter as a function of the radial distance from the ionizing stars. These re- 
sults show the reliability in using the R23 index to characterize the metallicity of H II regions 
even when only a fraction of the total area is covered by the observations. 
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1 INTRODUCTION 

Emission-line spectra of Galactic and extragalactic HII regions 
are normally obtained to characterize their physical properties, 
to extract information of the stellar populations that ionize them 
dOsterbrock & Ferlanoj [2006h as well as to study the variation 
of the physical p roperties of the star-forming regions across the 
galaxy disc (e.g. iMcCall. Rvbslri & Shields! 1 19851 ; IVflchez et ail 
ll988l : IVila-Costas & Edmunds! 19921) . The H II regions are far from 
the idealized spherical ionized gas clouds with homogenous phys- 
ical properties, as has clearly been shown in detailed studies of 
the nea rest Galactic and ex t ragalactic Hll regions: e. g. Orion 



Nebu la dBaldwin etal . 1991; Kennicutt et al. 2000; Sanche z et ail 
l2007h . 30 Doradus dMathis et al.ll 19851: iKennicutt et alj|2000h and 
NGC 604 dGonzalez-Delgado & Perez!l2000l : lMaiz-Apellaniz et al.l 
120041) . among others. The derived physical properties, generally 
obtained from long-slit spectroscopic measurements centred on 
the most intense knots, are not necessarily representative of the 
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condi t ions in all the l ocations within the regions dOev & Shields! 
l2000l : IOevetai]|2000h . The variations of the physical properties 
are especially seen in star-forming regions where the action of 
the massive stars strongly in fluence the surrounding medium (e.g. 
Kobulnickv & SMlmar] [l996l . 1997). Moreover, in some cases the 
gas within the H II regions is ionized by multiple star clusters which 
are not always centrally located. This will produce different tem- 
perature and ionization structures t han the normally assumed con- 
figurations from spherical mo dels dErcolano. Bastian & Stasinskal 
l2007l ; |jamet & Morissej|2008l) . 

In order to study the variations of the physical properties 
within the H II regions, the usual technique carried out up to now 
consists of a set of spectroscopic observations with slits located 
at different positions in the region and covering as much surface 
as possible. This technique has limitations however: it requires a 
large amount of observation time and does not cover, in general, 
the complete face of the region; thus interpolations of the observa- 
tional parameters in the gap between slits need to be made. Integral 
Field Spectroscopy (IFS) overcomes these limitations and offers 
the opportunity of extracting spectroscopic information at different 
positions across a continuous field of view. Therefore, this kind of 
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surveys allows the study of the variation of the physical properties 
within star-forming regions. 

We present here IFS observations of NGC 595, the second 
most luminous Hll region in M33. NGC 595 presents an angu- 
lar size of ~1 arcmin, which at t he distance of M33 (840 kpc; 
iFreed man. Wilson&Madorelll99lh translates into a linear physi- 
cal size of ~200 pc. This makes the region particularly suitable for 
being studied with the IFS technique. Although the current IFS fa- 
cilities usually have a relatively small field of view, the proximity 
of this region (which implies a relatively high surface brightness) 
makes possible to map it in a reasonable amount of observing time. 

NGC 595 has an Ha shell morphology that shows the ac- 
tion of the stellar winds of the massive stars located in its interior . 
The stellar content studied by iMalumuth. Waller & Parked J 19961) 
using optical photometry reveals the existence of ~250 OB-type 
stars, ~13 supergiants and 10 candidate Wolf-Rayet (WR) stars. 
These authors derived an age of 4.5 ±1.0 Myr for the stellar 
cluster, which is consistent with the age predicted from the syn- 
thesis of integrated sp ectra in th e far-ultraviolet (FUV ) wavelength 
range dPellerirj [2006). Recently, iDrissen et al.l d2008l) spectroscop- 
ic ally confirmed nine of the WR candidates previously identified 
bv lDrissen, Moffat & Sharaldl993l) using photometric observations 
and discarded two possible WR candidates from their sample. The 
physical properties for NGC 595 have been derived with long- 
slit spectroscopic observations located at the most intense knots 
within the region (Vilchez et al. 1988 ) and recently w ith echelle 
spectroscopy dEsteban et alj20091) . IVflchez et alj ( 1988) obtained a 
temperature of T c ~ 8000 K, an electron density consistent with 
the low density limit and a metallicity of 12 + log[ Q/H] = 8.44 ± 
0.09 (Z = O.6Ze0), which was recently confirmed by Este ban et al.l 
(2009). The extinction within the region has been stu died previ- 
ously by other authors: iBosch. Terlevich & Terlevicr] d2002l) ob- 
tained the extinction using the Ha/H/3 emission-line ratio and 
Viallefond, Donas & Gossl dl983l) derived the extinction with the 
Ha and thermal radio emission ratio. Relano & Kennicut tl J2009I) 
revised the extinctions using radio-to-Ha emission and derived new 
values using the 24 /im-to-Ha integrated emission ratio, finding 
consistent results from both methods. They also found differences 
in the emission distributions at several wavelength ranges: while the 
infrared emission at 24 /im and Ha are spatially correlated with 
each other, the FUV emission is located within the observed Ha 
shell structure of the region. They suggest that the dust emitting at 
24 fim will probably be mixed with the ionized gas of the region 
and will be heated by the central ionizing stars. 

In this paper, we analyse IFS observations covering the whole 
face of NGC 595. We are able to perform a detailed analysis of the 
principal emission-line ratios which are widely used to describe 
the main properties of H II regions such as extinction and density 
structure, ionization parameter, metallicity, existence of shocks and 
evolutionary state. With these observations, we are in a position to 
test the influence of the geometry and the stellar distribution within 
the region on the integrated fluxes of the emission lines, as has been 
predicted by the models. The comparison of our results with those 
obtained from long-slit spectroscopy makes it possible to estimate 
the bias of the observations when only one part of the region, gen- 
erally the most intense, is covered. Finally, the power of the IFS 
and the quality of the observations presented here allow us to anal- 



1 We assume solar abundance from Asplund, Grevesse & Sauval ( 2005), 
12 + log[Q/H] = 8.66 



yse the WR stellar content of the total surface of the region and to 
make a complete census of this stellar population in NGC 595. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations 

The data were obtaine d with the Potsd am Multi Aperture Spec- 
trophotometer, PMAS dRofh et al]| 20051) at the 3.5 m telescope in 
Calar Alto on the night of 2007 October 10. We used the lens array 
(LARR), made of 16 x 16 2 elements, with the magnification scale 
of 1.0 x 1.0 arcsec 2 . This configuration allows o bservation at hig h 
spatial resolution and with a filling factor of 1 .00 ( Kelz et al. 2006). 
We used the V300 grating, which provides with relatively low spec- 
tral resolution (3.40 A pix -1 ) but allows to cover the whole optical 
spectral range (3650-6990 A). 

Given the relative small field of view of the LARR in com- 
parison to the angular size of NGC 595, we made a mosaic of 13 
tiles to map the whole H II region, covering a field of view of 47 x 
92 arcsec 2 , corresponding ~174 x 340 parsec 2 at the distance of 
M33. The distribution of the different tiles is shown in Fig.[T]over- 
plot on an Ha emission-line image from N OAO (Nationa l Optical 
Astronomy Observatory) Science Archive dMassev et alj|2007l) . In 
order to have all the data in common units, we made a 2.0 arcsec 
overlapping between contiguous tiles. 

We obtained 3-4 exposures of 400 s each per tile, depending 
on the transparency of the sky and the surface brightness of the 
mapped region. Under non-photometric conditions, seeing ranged 
typically between 1.2 and 1.8 arcsec, while the air masses ranged 
between 1.9 and 1.0. In addition to the science frames, continuum 
and HgNe arc lamps exposures as well as sky background frames 
were obtained interleaved among the science frames. This allowed 
us to minimize the effects due to flexures in the instrument and 
ensure a proper sampling of background variations. 

Finally, exposures of the spectrophotometric standard-stars 
BD+28D4211, HR153, G191-B2B were obtained in order to cor- 
rect for the instrument response and perform a relative flux calibra- 
tion. 

2.2 Data reduction 

The data reductio n fo llows the process described in 
lAlonso-Herrero et al.l d2009h with some minor modifications. 
We performed the basic reduction of each individual frame using a 
set of homemade scripts created under the IRAF enviromenQ The 
PMAS's CCD presents some structure, so bias was removed using 
a masterbias made out the combination of 10 bias frames. We then 
identified the location of the spectra in the CCD using continuum 
lamp exposures (i.e. tracing) and extracted them. 

We performed a wavelength calibration using the HgNe lamp 
exposures. In order to estimate the quality of this calibration, we 
fitted a Gaussian profile to the [Ol]A5577 skyline in each spec- 
trum. The standard deviation for the centroid of the line in each 
individual exposure ranged typically between 0.1 and 0.4 A. We 
then corrected from the lens to lens sensitivity variations by creat- 
ing response images from the continuum lamp exposures. 

2 The Image Reduction and Analysis Facility IRAF is distributed by the 
National Optical Astronomy Observatories which is operated by the asso- 
ciation of Universities for Research in Astronomy, Inc. under cooperative 
agreement with the National Science Foundation. 
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Figure 1. Mosaic to map NGC 595 overplotted on a continuum-substracted 
Ho direct image from NOAO (Natio nal Optical Astronomy Observatory) 
Science Archive iMassey et al. 2007). The orientation is north up and east 
to the left. The Ha image is shown in logarithmic stretch to better enhance 
all the morphological features of the H II region. 



At this stage, cosmic rays were removed in each individual 
fram e using the spectro scopic version of the L. A. Cosmic rou- 
tine ( van Dokku mf l200lh . Input parameters were chosen in order 
to eliminate as most cosmic rays as possible. In the few possible 
cases where some spectral ranges close to emission lines were af- 
fected by cosmic rays (namely Ha and, in some particular cases, 
H/3 and [O III]A5007), these were removed manually. 

The next step was the relative flux calibration performed us- 
ing the G191-B2B standard star exposure, which was observed at 
very low airmass (i.e. 1.06). The flux calibration was done using 
the IRAF tasks standard, sensfunc and calibrate. In 
order to estimate the uncertainties associated to this calibration, we 
made a comparison of the sensitivity functions derived for G191- 
B2B and for the other two standards. For the spectral range of 

4500 7000 A, uncertainties were ~5 per cent, while in the 

bluer part of the spectra (i.e. <4500 A), they could reach 15 per 
cent. Note, however, that since it was not strictly necessary for the 
goals of the present work and the night was not photometric, we 
did not perform any absolute flux calibration. 

For sky subtraction, we created a high signal-to-noise (S/N) 
spectrum for each of the interleaved sky exposures by combining 
its 256 spectra. This collection of combined spectra sampled the 
sky variations along the night. The sky in each individual tile was 
subtracted using the sky spectrum obtained within ~ 1 h of the time 
of the object exposure. Because of the non-photometric conditions, 
small scaling factors (i.e. between 0.8 and 1.5) had to be included 
in the sky spectra in order to minimize the skyline residuals. 

Finally, we made a unique data cube for the whole Hll 
region by mosaicking the individual frames using the offsets 
commanded at the telescope. For this purpose, we u sed the 
create jnosaic . pro routine from the P3d package dBeckerl 



l2006t iRoth et al.ll2005h which corrects from the effect of the dif- 
ferential atmosph erical difraction by using the expresion given in 
Filippenkl] jl982h . All the frames were put to common units by 
using the overlapping bands between the individual tiles. 



2.3 Line fitting and map creation 

We used MPFITEXPR algorithrr0 implemented by C. B. Mark- 
wardt in IDL t o fit the emission lines in each spatial element (here- 
after, spaxel) (Markwardt 2009). The advantage of this algorithm 
is the simplicity of imposing restrictions in the parameters of the 
fit, which is carried out using Levenberg-Marquardt least-squares 
minimization. 

For each spaxel we selected the wavelength range that in- 
cludes the emission lines we are interested in fitting and then per- 
formed the minimization of the observed profile to a Gaussian 
profile describing the emission line plus a one-degree polynomial 
function for the continuum subtraction. As an illustration, for the 
Hq+[N II] emission lines we selected a wavelength range of 6500- 
6650 A, imposed a fixed wavelength separation between the Ha, 
[Nll]A6548 and [Nll]A6584 given by the redshift provided by 
NEE0 and performed the fit keeping the same linewidth for all 
the emission lines and a nitrogen line ratio of 3. We estimated the 
noise of the corresponding spectrum using the standard deviation 
of the adjacent continuum. The noise defined in this way was used 
to obtain the S/N for each emission line fitted. We performed the fit 
in each spaxel and rejected in the procedure those fits with a S/N 
smaller than 5. The remaining spectra were visually inspected and 
bad fits were interactively rejected. 

The same procedure was applied to all the observed emis- 
sion lines. The flux for each emission line was obtained by inte- 
grating the area defined by the best-fitting Gaussian profile. Then 
we used the derived flux and the position within the data cube for 
each spaxel to create an image (a map) that can be treated with 
standard astronomical software. We performed the astrometry of 
the emission-line maps using a continuum map close to the Ha 
emission line derived in the fitting procedu re and a R-Band im age 
obtained from the NOAO Science Archive (Massey et al. 2006) to 
identify the most intense stellar clusters in the continuum map. The 
accuracy of the astrometry is better than 1 arcsec, the angular size 
of the spaxel in our field of view. 



3 INTEGRATED PROPERTIES 

Before discussing the emission-line maps of NGC 595, we first 
analyse the spectrum for the whole Hll region. In Fig. [2] we show 
the integrated spectrum in arbitrary units obtained by coadding 
the signal from the ~3000 spaxels. We use two different nor- 
malizations in the figure to show all the observed emission lines. 
We checked our method (Section |23) by comparing the measured 
fluxes with those derived using the task splot from IRAF. We 
find differences in the fluxes of ~ 1 per cent, which shows that our 
method is reliable and gives consistent results. We corrected the 
fluxes using the observed Ha/H/3, Hy/Hff, HS/HQ, an d He/H fi ra- 
tios and the mean extinction law o f ISavage & Mafhisl jl979h . The 
reddening coefficient C(H/3) was derived in an iterative procedure 



3 See http://purl.com/net/mpfit 

4 http://nedwww.ipac.caltech.edu 
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Table 1, Dereddened emission-line ratios with respect to / (H/3) detected in 
the integrated spectrum for the complete field of view of our observations 
(Fig.EJ. The errors are a combination of the profile-fitting procedure and 
the uncertainty in the co ntinuum subtraction. f(A ) is the reddening function 
normalised to H/3 from Savage & Mathis] ll979l) . 



Emission Line 


/(A) 


/(A) 


3727 [O II] 


0.252 


2.86± 0.05 


3889 H8+Hel 


0.220 


0.12± 0.02 


3970 He + [Ne III] 


0.204 


0.11±0.02 


4102 H<5 


0.177 


Q.21± 0.01 


4471 Hel 


0.100 


0.04± 0.01 


4340 H7 


0.131 


0.44± 0.03 


4959 [Om] 


-0.030 


0.295± 0.006 


5007 [Om] 


-0.043 


0.95± 0.01 


5876 He I 


-0.220 


0.099± 0.003 


6563 Ha 


-0.314 


2.94± 0.02 


6584 [Nil] 


-0.319 


0.505± 0.005 


6678 Hel 


-0.329 


0.032± 0.003 


6717 [S 11] 


-0.335 


0.267± 0.005 


6731 [S 11] 


-0.336 


0.1 86± 0.004 



assuming the same EW a b s for all the Balmer li nes and the theoreti- 
cal ex pected Balmer ratios at T c = 7.6 x 10 3 K dStorev & Hummed 
1 19951) . Following this method, we find a coefficient of C(H/3) = 
0.17 ± 0.03 and estimate the stellar absorption equivalent width to 
be always EW a b s ^2 A. 

The dereddened emission line ratios with respect to /(H/3) 
are given in Table [7] The corresponding errors are a combination 
of those derived in our profil e-fitting procedure and the unce rtainty 
in the continuum subtraction (Gonzalez-Del gado et al.ll 994). They 
do not include the uncertainty in the calibration that can be between 
5 and 15 per cent, depending on the exact wavelength range (see 
Section 2.2). We find d ifferences of ^20 per cent with the emission- 
line ratios reported by Esteban et al. (2009), except for the follow- 
ing lines: [On]A3727, [Sn]A6717, [Sn]A6731 and H8+Hel for 
which we find differences of 40 per cent for [O II]A3727, 100 per 
cent differences for the [S II] emission lines and 24 per cent for the 
last case. The differences for [Oll]A3727 and H8+Hel can be at- 
tributed to the diffe rent spectral resolutions of both observations: 
Este ban et all (2009) are able to resolve the two [O II] lines and the 
H8 and Hel separately, while the fluxes reported here are the com- 
bination of the unblended emission lines. The differences in the 
[S II] emission lines are related to difference s in the apertu r e used 
in each study: the slit field of view used by lEsteban et alj (2009) 
(5.76 x 1.7 arcsec 2 ) is small and it is located over a knot with high 
Ha surface brightness, while the fluxes reported here correspond 
to the t otal area of the regio n. Integrating over the slit field of view 
used in lEsteban et ail 12009). we obtain differences of ~ 10 per cent 
for both [S II] emission lines. As we will see in Section 14.41 the 
[S II] emission-line flux in each spaxel shows a strong dependence 
on the same spaxel Ha flux, whic h explains the differ ences found 
here with the [S II] fluxes given in lEsteban et all J2009h . 

In Table [2] we show the principal diagnostic emission-line 
ratios derived from the integrated spectrum as well as the phys- 
ical parameters derived from them. The integrated reddening co- 
efficient corresponds to an extinc tion of 0.36 mag at Ha, w hich 
agrees with the values reported in lRelafio & Kennicuttl (2009) de- 
rived from other methods. 

Due to the contribution of a strong sky HgI4358 emission line 



and the low intensity of [Om]A4363 we were not able to detect 
this emission line, and thus a measurement of the temperature with 
these observatio ns cannot be mad e. We report here the mean value 
obtained from lEsteban et al.l 12009J) using high-resolution echelle 
spectroscopy, which has been used to obtain the electron density 
from the [S II]A6717/[S II]A6731 emission line ratio. We derived 
a nominal density value of 13 cm -3 , for a line ratio of 1.43 (see 
Table O, whi ch agr ees with the low density scenario derived by 
lEsteban et al.l j2009h . Taking into acco unt the extinction-corrected 
Ha luminosity of NGC 595 reported in Relano & KennicutJ b009T) 
and the electron densities derived here, we estimate an ionization 
parameter of q c ff=1.9xl0 7 cm s _1 for an electron density face 
value of n c =20 cm -3 . 

The metallicity of the H II region can be derived using infor- 
mation from several emission-line ratios. For this purp ose, we will 
use a s metallicity calibrator the widely used R23 index I Pag el et all 
1979). As we will discuss later in Section 1431 this parameter is 
substantially constant across the face of the region, and thus we use 
it here to derive a representative value for the Hll region metal- 
licity. The calibration of R23 with metallicity has been extensively 
studied in the literature, either using p hotoionization models (e.g. 
lMcGaughTl 991: Kewley & Dopita 2002) or from empirical calibra- 
tions (e.g. |Pilvuginl l2000: Pilyugin 2001). One of the major prob- 
lems of R23 is that it is double valued; thus an independent metallic- 
ity diagnostic must be used in order to set the region in the upper or 
lower branch in the Z-R23 diagram. We use the observed [N Il]/Ha 
line ratio to make a first estimation of Z : follo wing the relation 
given by iDenicolo, Terlevich & Terlevichl (120021) we obtain 12 + 
log[0/H] = 8.56, while the calibration of fpettini & Pagell J2004l) 
gives 12+log[0/H]=8.41. Both results shows that we are able to 
use the upper branch of the Z -R23 diagram. We apply two cali- 
br ations to obtain Z from R23, McGaughl <fl99 lh as parametrised 
by Kobulnicky, Kennicutt & Pizagno (1999), which relie s on pho- 
toioni zation models, and the empirical calibration of IPilvuginl 
(2001). The first one gives 12 + log[0/H] = 8.75, while the em- 
pirical relation yields 12+log[0/H]=8.43. The mean value between 
these two calibrations, 12+log[0/H]=8.59, was finally chosen to 
characterize the oxygen abundance for the whole H II region, with 
a typical uncertainty of ±0.32 dex taking into account the different 
results from both calibrations. This value agrees with the range of 
previous values given in th e literature: 8.45 ± 0. 03 (t 2 = 0.000) and 
8.69 ± 0.05 (t 2 = 0.0 36) bv lEsteban et"ail J2009h and 8.44±0.09 by 
IVflchezetallh988l) . 



4 EMISSION LINE MAPS 

4.1 Structure of the ionized gas and the stellar component 

In Fig. [3] we show a map of the observed Ha flux for NGC 595 
derived with the procedure explained in Section 12.31 The simi- 
larity with the observed Ha flux obtained from direct image (see 
Fig. HJ shows the power of the IFS observations to reproduce the 
surface brightness distribution for this particular region. The Ha 
shell structure is clearly delineated, it extends from the south-west 
to the north in an arched structure with Ha maxima located close to 
the ionizing stars. The separation between the central star clusters 
and the Ha maxima is ~6 arcsec, ~22 pc. The external low Ha 
surface brightness shell, seen in Fig. [T] eastwards to the brightest 
one, is also observed in the Ha map, which shows the high qual- 
ity of the data to trace the diffuse Ha emission of the region. In 
Fig. [3] we overplot emission contours of the Wide Field Planetary 
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Figure 2. Integrated spectrum of NGC 595 obtained by coadding the signal of all the spaxels in the field of view with two different normalizations to better 
visualize the different observed emission lines. Fluxes are in arbitrary units, and the main detected emission lines are listed in Table [T] Some residuals at 
~559() A and ~63()0 A produced in the sky subtraction procedure under non-photometric conditions are seen in the integrated spectrum. 



Table 2. Principal diagnostic emission-line ratios and integrated physical 
properties of NGC 595. Except for the electronic temperature that was taken 
as a mean value of the temperatures derived by Esteban et ah ( 2009J), the 
rest of the physical properties have been estimated here using our integrated 
fluxes (see te xt for details ) R 2 3= ([O II] A3727 + [O III] AA4959,5007)/H,3, 
as defined bv lPagel et all ll979l) . 



Parameter 


Value 


log [Nll]6584/Ha 


-0.765 ±0.005 


log [S II] (6717+673 l)/Ha 


-0.813±O.0O7 


log [O III] A(4959+5007)/[O II] A3727 


-0.361 ±0.009 


log R23 


0.614±0.006 


[Sn]A6717/6731 


1.43±0.04 


C(H/3) 


0.17±0.03 


n c (cm -3 ) 


<220 


q c ff (cms" 1 ) 


1.9xl0 7 


T c (K) 


7670± 116 (adopted) 


12+log[Q/H] 


8.59±0.32 



Camera 2 (WFPC2)/F336W image from the Hubble Space Tele- 
scope (HST) Multimission Archive at the Space Telescope Science 
Institute (MAST) Archive. The knots correspond to the location of 
the stellar clusters; two of the most intense ones are located in the 
centre of the shell and do not correspond to regions of high Ha 
emission; the third intense knot is within the Ha shell structure, 
slightly north-east from the others. Farther away from the most in- 
tense knots, there are others spatially correlated with diffuse Ha 
emission. 



4.2 Extinction in NGC 595 

We used the observed Ha/H/3 ratio to derive a map of the red- 
dening coefficient C(H/3) and to correct the emission-line fluxes 
at each spaxel of the field of view. The ratio was compared with 
the theoretical expected value of the Balmer dec rement for case B 
of rec ombination theory at T c = 7.6 x 10 3 K ( Storey & Hummei 
1995), and then by using the extinction law of ISavage & Mathis 
( 1979) we obtained C(H/3). The Balmer emission lines can be af- 
fected by the absorption of the underlying stellar component, which 
would affect the derived values of the reddening coefficient. We ex- 
plore the influence of this effect in our data by visual inspection of 
the spectra and by a comparison of the observed H/3 equivalent 
width (EW(H/3)) in emission and the expected absorption equiva- 
lent width (see Section^. Except for a circular region of ~3 pixels 
associated with the location of the central clusters where we mea- 
sure EW(H/3) of ~5-10 A, we find EW(H/3) within a range of 40- 
450 A. The effect of the underlying stellar population is negligible 
at the locations of the high equivalent widths. A closer inspection 
of the Balmer emission-line profiles at the location of the central 
clusters would suggest the presence of absorption signatures for 
H<5 and He which are marginally detected (S/N^3) and consistent 
withEWabs ~( 1.0- 1.4) A. 

A C(H/3) map is shown in the upper panels in Fig. [4] It 
presents a very concentrated distributi on within the region, s imilar 
to the extinction structure obtained bv lViallefond et alfdl983h from 
radio and Ha observations. We derive values for the extinction at 
Ha in the centre of the region within the range of 1.0-1.5 m ag, 
which agrees with the range reported in Viallefo nd et alj d 1983b for 
the Balmer extinction in the core of the nebula, (0.88-1.15) mag. 
With the reddening coefficient map we can then correct the ob- 
served spectrum in each spaxel within the region and obtain maps 
of the emission-line ratios that can be used to study the physical 



6 M. Relano, A. Monreal-Ibero, J. M. Vilchez and R. C. Kennicutt 



C(H/3) C(H/5) 




R.A. (2000) R.A. (2000) 



Figure 4. Upper: Map of the reddening coefficient C(H/3) for NGC 595 with 24 fim (left) and 8 (im (right) emission contours overplotted. The intensity 
contours for the 24 and 8 fim emission are at (2, 5, 10, 20, 40, 60, 80, 95) per cent of the maximum intensity within the region. A 1 per cent contour level 
corresponds to 3cr and ltr for the 24 and 8 fim emissions, respectively. Lower: Absorbed Ha luminosity of NGC 595 with 24 fim (left) and 8 fim (right) 
emission contours overplotted. The contour levels are the same as in the upper figures. 
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R.A. (2000) 

Figure 3. Map of the observed Ha flux derived using the fitting procedure 
described in Section 2.3. The field of view of the observations is depicted in 
Fig-H Each spaxel has a 1 X 1 arcsec 2 size and fluxes are in arbitrary units. 
The contours correspond to the WFPC2/F336W filter image obtained from 
the HST MAST Archive. The most intensity contours show the location of 
the ionizing stellar clusters within the region. 



properties of the Hll region. As an illustration, we show in Fig. [6] 
some of the emission-line ratio maps that will be analysed later in 
this paper. 

In Fig. [4] (upper panels ) we overplot contours of the 24 and 
8 /im emissions from Spitzer ( Relano & Kennicuttl2009l) . The loca- 
tion where C(H/3) is high corresponds within the spatial resolution 
(~1.5 arcsec for our observations and ~6 arcsec for the observa- 
tions at 24 ^m) to high values of the 24 /im emission. The emission 
at 8 fim also correlates to the maximum values of C(H/3), but it is 
not as concentrated as the emission at 24 fim (Fig. [4] upper left). In 
the lower panels of Fig. [4] we show the relation of the emission at 
24 and 8 fim with the absorbed Ha luminosity of the region. The 
Ha luminosity absorbed in the region is obtained as the difference 
between the total extinction-corrected Ha luminosity and the Ha 
luminosity corrected for the foreground Galac tic extinction (E(B- 
V)=0.041 dSchlegel. Finkbeiner & Davislll998l) is used to account 
for the Galactic extinction). As it is shown in Fig. [4] (lower left 
panel), the 24 fim emission spatially correlates with the absorbed 
Ha luminosity, with the maximum at 24 /im coinciding with the 
central position of the absorbed Ha structure. The spatial correla- 
tion of the 8 fim emission and the absorbed Ha luminosity is not so 
strong; the maximum at 8 fim is shifted outwards of the absorbed 
Ha shell structure. 

These spatial correlations are better shown in Fig. [5] where 
we plot the emissions integrated in elliptical concentric annuli as 
a function of the radial distance measured along the major axis of 
an ellipse centred at the location of the main stellar clusters (we 
take this position to be 2 arcsec westwards from the most intense 



stellar cluster marked as a red cross in the emission-line maps of 
Fig.(6]l. The major to minor axis ratio of the ellipse i s derived using 
the inc lination angle of the galaxy (i=56° for M33 Ivan den Berghl 
(2000)), and we choose a position angle of 129° for the major axis 
since this orientation better traces the shell structure of the region. 
In this configuration, we use rings of 2 arcsec widths to obtained 
the elliptical profiles. 

The Ha/H/3 emission line ratio was derived for each annulus 
in the original data cube. The resulting spectrum for each annulus 
was fitted in the same way as we explained in Section 12.31 and 
the corresponding fluxes at Ha, H/3 were used to obtain the radial 
profiles shown in Fig. a )- The absorbed Ha luminosity profile 
(Fig- SO was obtained in the same way we derived the absorbed 
Ha luminosity map. 

The comparison of Fig.[5]a) and (c) shows that the radial dis- 
tribution of the Ha/H/3 ratio and the 24 fim emission are similar, 
with the same radial decline towards larger radial distances and 
maxima slightly shifted (~2 arcsec, lower than the resolution of 
the 24 /im observations (~6 arcsec)). The comparison of Figs. Ob) 
and (c) shows that the distribution of the absorbed Ha luminos- 
ity and the 24 fim emission are also similar, with their respective 
maxima located at the same radial distances. The spatial relations 
shown here between the 24 fim emission, the extinction suffered by 
the ionized gas and the absorbed Ha luminosity support the idea 
that the 24 fim emission is produced by heating of dust probably 
mixed with t he ionized gas inside t he re gion, which was previously 
suggested by Rel ano & Kennicutn (2009). 

The emission at 8 fim also correlates to the maximum values 
of C(H/3), but the spatial distribution of the 8 fim emission does 
not follow the C(H/3) structure as the 24 fim emission distribution 
is (Fig. [4] upper panels). In Fig.[5ld), we plot the 8 fim emission 
elliptical profiles. The maximum at 8 fim is ~4-5 arcsec shifted 
towards higher radial distances than the maximum of the Ha/H/3 
ratio. This shows that the dust responsible for the Balmer extinc- 
tion is the dust emitting at 24 fim rather than the dust emitting at 
8 fim. Relano & Kennicutt (2009) showed that the 8 fim emission 
for NGC 595 is more related to the location of the main CO molecu- 
lar clouds identified in this region than to the position of the ionized 
gas. 

4.3 Density Structure 

We study the electron density using the [S II]A6717/[S II]A6731 
emission-line ratio. A map of this ratio is shown in Fig. [6] 
(bottom-right). The emission-line ratio ranges from 1.2 to 1.8 
(n c < 220 cm -3 ), showing that the Hll re gion has in general a 
low electron dens i ty, as was also found by lEsteban et a 1. (2009) 
and IVflchez et alj (l988h . The map shows no particular density 
structure, and further maps obtained from binning the data cube 
in 2 x 2 and 4x4 spaxels do not show any density structure neither. 
lLagrois & Joncasl j2009h show a density map for NGC 595 which 
extends in a wider field of view than ours; the lower density values 
in their map are associated to the location of HI neutral gas, while 
within the area cover by Ha emission their map shows no strong 
variations with density values of ~60-150cm -3 . 

However, the pronounced Ha shell morphology of the neb- 
ula would suggest a density enhancement at the location of the Ha 
maxima. Density variations within the H II regions have been re- 
ported before and in some cases, a r elation between the density and 
the H a emission has been observed I Casta neda. Vflchez & Copettil 
1992). In order to check carefully any possible density variation 
with the Ha emission, we have compared the radial profiles of 
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Figure 6. [S II]A6717,31/Ha (upper left), [N II]A6584/Ho (upper right), [O III]A5007/H/3 (bottom left) and [S II]A6717/[S II]A6731 (bottom right) emission- 
line ratio maps for the whole face of NGC 595. Only spaxels with emission-line ratios having relative errors < 30 per cent are shown. Extinction correction 
for the fitted emission lines was performed in each spaxel prior obtaining the emission-line ratios shown here. The red cross marks the location of the central 
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Figure 5. Radial profiles (from top to bottom) of the (a) Ha/H/3 emission 
line ratio, (b) absorbed Ha luminosity and Spitzer infrared bands at (c) 
24 pm and (d) 8 pin. For panels (a) and (b), we isolated the signal coming 
from concentric elliptical rings of 2 arcsec width projected onto the mayor 
axis of the ellipse and derived the corresponding integrated spectrum at each 
annulus. The center of the ellipse (-2,0) is expressed in the units of the 
emission-line maps shown in Fig. [6] We fitted the spectra in the same way as 
explained in Section l2~3l and derived the integrated emission-line fluxes for 
each concentric annulus (see text). For the 24 and 8 fim emission profiles, 
we plot the integrated fluxes for each elliptical annulus. The dot-dashed line 
represents the location of the maximum of the absorbed Ha luminosity. 
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Figure 7. Radial profiles of the observed Ha luminosity (a), 
[S II]A6717/[S II] A6731 (b), obtained in the same way as the profiles 
shown in Fig. [5] 



the [S II]A6717/[S II]A6731 line ratio (Fig.[7]b) with profiles of ob- 
served Ha emission (Fig.[7]a). 

The mean value for the [S II]A6717/[S II]A6731 line ratio in 
all the rings is (1.43±0.06), corresponding to the low density limit. 
The elliptical profiles show that there is no electron density varia- 
tion with the Ha luminosity distribution for this region. In this sit- 
uation, at the location of the maximum in Ha luminosity, which 
would correspond to an enhancement of mean electron density, 
there should be an increase of the filling factor producing a higher 
emission measure along the line of sight. In order to further study 
this effect, detailed photoionization models of the region would 
need to be performed (Perez-Montero et al., in preparation). 

4.4 ionization structure 

Emission line ratio maps of [S II]A6717,31/Ha, [Nll]A6584/Ha, 
[Om]A5007/H£ and [Sn]A6717/[Sn]A6731 are shown in Fig.© 
The last ratio gives information of the density of the H II region and 
was studied in previous section; the others allow us to study the ion- 
ization structure of the region. Lower values of [S II]A6717,31/Ha 
and [N II]A6584/Ha are located close to the position of the ion- 
izing stars, corresponding to the high excitation zone within the 
region, while at the outskirts the values of these ratios are higher, 
depicting the low excitation zone. The same effect is seen in the 
[O IIl]A5007/H/3 map; the high excitation zone closer to the central 
stars is shown by the high values of this ratio, while the low excita- 
tion zone is observed far from the central region, corresponding to 
low values of [O lii]A5007/H/3. 
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Using BPT diagrams ( iBaldwin. Phillips & Terlevicbl 1 198 ll ; 
IVeilleux & O sterbrock Il987h and the emission-line ratio maps in 
Fig. [6] we can study whether non-photoionizing mechanisms (e.g. 
shocks) are important within NGC 595. In Fig. [8] we show the 
BPT diagrams for all the spaxels in the field of view with mea- 
surements of [O UiyU/3, [N Il]/Ha and [S Il]/Ha grouped into bins 
of Ha flux. We also show the separation between active galac- 



different studies i 


Veilleux & Osterbrocklll987l;lKewlev et alj:2001; 


Kauffmann et al. 


2003; Stasinska et al. 200rj). As can be seen in 



the figure, all the spaxels lie below the separation lines, which 
shows that the main mechanism producing the observed emission 
lines within NGC 595 is probably photoionization. We have com- 
pared (not shown here) the location of the spaxels in these diagrams 
with the one for a wide sample of H II regions from the literature. 
The area occupied by the set of Hll regions from the literature 
could be considered as the representative location for H II regions, 
where photoionization dominates. The comparison shows that all 
our spaxels fall into this area; therefore, if non-photoionizing mech- 
anisms are present in the region, these are not contributing substan- 
tially to the emission of NGC 595. However, a detailed comparison 
of these obse rved line rat i os for some spaxels with the predictions 
of models by I Allen et al] (2008) suggests that somewhat localized 
effects of shocks can not be disregarded. 

In the top panels of Fig. [8] we study the [O IIl]/H/3 ver- 
sus [N Il]/Ha, while the bottom panels show [O IIl]/H/3 versus 
[S Il]/Ha. The distribution of the data is different in each set of 
panels. In the case of [O IIl]/H/3 versus [N Il]/Ha, the spaxels with 
a high Ha flux (top-left hand panel) correspond to those showing 
high excitations levels (high values of [O IIl]/H/3); however, the low 
Ha flux spaxels (top-right hand panel) cover the whole range of ex- 
citations within the nebula. The highest [O IIl]/H/3 values in the low 
Ha flux distribution are in the x-axis range of log([N Il]/Ha)^ (- 
1.3,-1.0). These spaxels are located close to the stars (north-west 
of the central clusters), where no high intensity Ha emission is ob- 
served (see top-right hand panels in Figs|3]and |6). 

The data distribution in the [O m]/H/3 versus [S Il]/Ha dia- 
gram (bottom panels of Fig. [8} shows a similar but no so strong 
trend as in the top panel. The separation between the distribution of 
the high Ha flux spaxels (bottom-left hand panel in Fig. [8} and the 
integrated value for the region, marked in the same plot as a black 
cross, is remarkable. This shows the effect of the aperture selection 
in the estimation of the [S II] fluxes and reflects the differences we 
have found between the result s for the integrated m easurements re- 
ported here and those given in lEsteban et alU 2009). only covering 
one of the most intense knots in the region (see Section[3]l. 

In Fig. [9] we compare t he observed BPT diagrams with mod- 
els from Dopita et al. (2006), which take into account the effect of 
the stellar winds on the dynamical evolution of the region. The ion- 
ization parameter is replaced by a new variable R that depends on 
the mass of the central cluster and the pressure of the interstellar 
medium ( R = (Mci/M )/(P o /fe), with P a /k measured in (cm -3 
K)). Using an estimate for the mass of the cluster of 2x 10 J Mq 
jRelano & Kennicutfc 009) and the temperature listed in Table|2] a 
density value of 13 cm -3 (Section 3) gives log R=0.3. The mod- 
els shown in Fig. [9] are for log R=0.0. All the spaxels fall in the 
area delimited by the lines at Z=Zq and Z=O.4Z0, which agree 
with the estimate of the metallicity of the region from other meth- 
ods (see Section[3j. The age of the region derived from the models 
(~3.0-3.5 Myr ) also agrees with previo us estimates in the l iterature 
(4.5± 1.0 Mvr iMalumuth et all Jl99^) and 3.5±0.5 Myr IPellerinl 
(2006)). The location of the points in the right-hand panel in Fig. [9] 



depends on the Ha flux: spaxels with high Ha flux have generally 
lower values of [S Il]/Ha, while those with low Ha flux are located 
towards the right in the x-axis. 

4.5 WR stellar population 

WR stars are very bright objects which show strong broad emis- 
sion lines in their spectra. They can be classified as nitrogen (WN) 
stars (those with strong lines of helium and nitrogen) and carbon 
(WC) stars (those with strong lines of helium, carbon and oxygen). 
They are understood as the result of the evolution of massive O 
stars which lose a significant amount of their mass via stellar winds 
showing the products of the CNO b urning first ( WN stars) and the 
He burning afterwards (WC stars) dContilll976h . The short phase 
for the WR stars makes their detection a very precise method for 
estimating the age of a given stellar population. Typically, an in- 
stantaneous burst of star formation shows these features at ages of 
~ 2 — 6 Myr and even at a m ore limited range at very low metal - 
licities dLeitherer et ai] Fl999). The presence of WR stars can be 
recognised via the WR bumps around A4650 A (i.e. the blue bump, 
characteristic of WN stars) and A5808 A (i.e. the red bump, char- 
acteristic of WC stars). 

The WR popul ation in NGC 595 has alre ady been studied in 
iDrissen et all ( 1 19931) and iDrissen et ail ( 120081) by means of high 
spatial resolution narrow imaging of the central area of the neb- 
ula, first, and spectroscopic follow-up, afterwards. They detected 
11 WR candidates: nine of them were spectroscopically classified, 
while those named as WR10 and WR1 1 were removed from the 
candidate list . The WR survey done by IDrissen et al. I dl993l) and 
IDrissen et alj d2008h only covers the central part of the region with 
a field of view of 30" x 35", and the spectral resolution of their ob- 
servations is ~9 A. The observations presented here have higher 
spectral resolution (3.4 A) and cover an four times greater. This al- 
lows us to make a complete census of the WR population of the 
region. 

In order to make a first identification of the WR stars, we 
constructed continuum-substracted images at 4490-4540 A and 
5645-5695 A spectral ranges. Fig. 1 10 ) s hows the positions of the 
WR features detected by Driss en et al.l d 19931) overplotted on our 
continuum-substracted 4490-4540 A map. There is a good corre- 
spondence within the errors of the astrometry. All the candidates 
but WR5 (which is not covered in the current data, as shown in 
Fig. [lOt and WR7 (the one with the faintest measured bump) are 
detected by visual inspection. In addition, we have detected a new 
candidate towards the north of the H II region. WR3 and WR1 were 
also detected in the map for the red bump (not shown). This map 
was less reliable since important residuals in the 5577 A sky line 
do not allow to select a reliable spectral range for the continuum 
subtraction. 

We extracted the integrated spectra of the WR candidates by 
coadding typically ~ 3 — 10 spectra centred at the position of the 
WR detections. Fig. [TJJ shows the spectra for our new candidate 
(botto m panel), as well as for the WR catalogued bv lDrissen et al.l 
(2008) (top panel), and illustrates the power of the IFS data to iden- 
tify and characterize WR stars. In all the spectra shown in Fig. 1111 
the blue bump characteristic of WN populati ons is clearly visible . 
We compare our spectra with those shown in lDrissen et al.1 d2008l) . 
These authors present the spectra of four WR stars, except for WR5 
(not covered by our observations) the spectra of the other stars 
( WR1, WR3 and WR9) show the same features as those presented 
bv lDrissen et alj d2008l) . We confirm the classification of the WR 
stars catalogued by these authors, all of them being WN stars ex- 
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Figure 8. BPT emission-line ratio diagnostics I Baldwin et al . 198h lVeilleux & Oster brock 1987) for all the spaxels in the field of view where we measure 
[O III]/H/3, [N II]/Ha and [S II]/Ha in bins of different Ha flux. The black cross is the value obtained using the total integrated flu xes of the H II region. The 
lines show the separ ation between A GN and NSF galaxies from dif ferent studies: continu ous line - fVeilleux & Osterbrock 1987), dot line - I Kewle v et alj 
1200 ll) . dashed line - fcauffmann et alj2003l) . and dashed-dot line - iStasiriska et al. 2 0061) . 
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Figure 9. BPT emission-line ratio diagnostics as shown in Fig. [8] compared with models from Dopita et al. (2006) (log [(M<-ji/MQ)/(P o /fc)]=0.0). In the 
left-hand panel we plot all the spaxels without making any distinction in flux, while in the right-hand panel we code the points for different bins of Ha flux. 
We plot the models for two different metallicities (Z=Zq and Z=0.4Zq) covering the range of metallicities reported in Section|3] We also mark the temporal 
evolution of the region at 2.5, 3.0 and 3.5 Myr. 
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Figure 10. 4490-4540 A map derived from the PMAS data with contours 
of the Ho emission (see Fig. [3J in logarithmic scale. Continuum has been 
substracted by averaging the spectral ranges of 4650-4750 A and 4755- 
4805 A. We have ov erplotted the positions of the WR candidates listed by 
iDrissen et alj [ 19931) with dark grey circles. Contours are in a logarithmic 
scale and over an interval of 1 .45 dex in steps of 0.36 dex. The current data 
show a new candidate towards the northern part of the nebula which has 
been indicated with a black circle. Orientation is the same as in the rest of 
the maps and fluxes are in arbitrary units. 



cept WR3, whose red bump is clearly observed at 5800 A, showing 
that it is a WC star. 

We estimate the EW for the blue bump as a sum of the indi- 
vidual EW for the emission lines seen in the 4600-4700 A wave- 
length range. Except for WR3 and WR4, the EW(blue bump) of 
the stars i n Fig. | ll | ranges w i thin 1 5-40 A, which agrees with the 
models o f Schaerer &Vaccal (l998) (see their fig. 1 1 for metallic- 
ities of Z=0.08 and Z=0.02). We derive an EWjblue bump)= 28 A 
for WR9, similar to the value reported bv lDrissen et al.l ( 120081) . For 
WR3 and WR4 we find EW(blue bump)^60 A, which are consis- 
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Figure 11. Extracted spectrum for the new WR candidate (bottom) and the 
WR stars catalogued by Drissen et al 1 fcoosh (top). Vertical arrows mark 
the position of the expected blue and red bumps. 



tent w ith p revious values given in the li terature jMassev & Contil 
dl983l) and lArmandroff & Massed dl99lh . respectively). 

The spectrum of our WR candidate is shown in the bottom 
panel of Fig. [TT] The blue bump is clearly seen in the spectra 
but no red bump is visible, which shows that this star is proba- 
bly a WN star. We derive an EW (blue bump)=\A A for our new 
candidate, which i s with in the expected range of values given in 
ISchaerer & Vaccal dl998h . We cannot rule out the possibility of 
this star being an Of rather than a WR star (R Crowther, pri- 
vate communication), but further spectroscopic observations at 
higher resolution are required to assess this issue. It is interest- 
ing to note the isolation of our new WR candidate at a distance 
of ~185 pc away from the location of the rest of the WR stars. 
Assuming a time of 3-4 Myr, consistent with the age of the H II 
region, we find radial velocities of 45-60 km s" 1 for the WR can- 
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didate. These velocities are in the range of v a lues expected for 
runaway WR and O stars dMoffat et all 1 19981 : iDrav et al.1 [20o3 ; 
iHoogerwerf. de Bruijne & de Zeeuwll2001 ). 



4.6 Metallicity, ionization parameter and structure of the 
region 

It is well known that the R23 index (R 2 3= ([O II] A3727 + [O III] 
AA4959,5007)/H/3) does depend on the ionization parameter but 
the question that arises is how strong is the variation of this 
parameter within NGC 595? Several observational studies have 
claimed evidence for a relative constant value of R23 in spite of 
the strong variations of the ionization parameter within H II re gions 
dOev et alj|200d. bey & Shieldsll2000l : Ikennicutt et aHl2000l) . Re- 
cently, lErcolano et al. U2007h suggest that the geometrical distribu- 
tion of the ionizing stars could account for the discrepancies in the 
metallicities derived using classical metallicity calibrators. In or- 
der to shed light on this problem, we have studied the variations 
across the surface of NGC 595 of several classical emission-line 
ratios proposed as metallicity calibrators and compared them with 
the variation of [O III] AA4959, 5007/[O II] A3727, a tracer of the 
ionization parameter. 

In the top panels of Fig. [T2] we show maps of the [O III] 
AA4959, 5007/[OII] A3727 (left) and R 23 (right). There is a clear 
radial trend of [O III] AA4959, 5007/[O II] A3727 with the radial 
distance from the central stars: close to the ionizing stars we find 
high values of this ratio, implying a high ionization parameter, 
while moving towards larger distances from the ionizing stars the 
emission-line ratio declines, showing that the ionization parameter 
is lower at these locations. A map of the R23 index for the whole 
region is shown in Fig. [12] (top right). We have quantified the vari- 
ations for both maps in the lower panels of Fig.Q/2] Using the same 
elliptical integration as we have explained in Section |4~2l we de- 
rived the values for [O III] AA4959, 5007/[O II] A3727 and R23 in 
elliptical annuli of 2 arcsec width. While there is a strong gradi- 
ent of 1 order of magnitude from the centre to the outer parts of 
the region for [Oin] AA4959, 5007/[On] A3727 (Fig. [12] bottom 
left), R2 3 (black diamonds in the bottom right-hand panel) shows 
a variation of 0.13 dex over the whole radial distance. This shows 
the reliability in using R23 as an emission-line ratio to obtain the 
metallicity of this H II region; even when only a part of the region 
is covered by the observations, which is normally the case for long- 
slit spectroscopy, the value derived for R2 3 can be used as a rep- 
resentative one for the Hll region. The variation of R23 in Fig. 1121 
(bottom right) would translate into a nominal difference of 0.17 dex 
in log [Q/H] in the case we would use the calibration o f lMcGaughl 
dl99ll) . 

For comparison, we have included in Fig. Q/2] (bottom right) 
other emission-line ratios proposed in the literature as metallic- 
ity calibrators: [Nil] A6584/[Oll] A3727 (blue diamonds), [Nil] 
A6584/[SII] AA6717, 6731 (green diamonds), [Nil] A6584/Ha 
(red diamonds), and [Nil] A6584/[0 III] A5007 (magenta dia- 
monds). The figure shows that [N Il]/Ha and [N Il]/[0 III] vary sig- 
nificantly with the radial distance. The variation is particularly sig- 
nificant for [N Il]/[0 III], 0.92 dex, which is exp ected since this ra- 
tio str ongly depends on the ionization parameter ( Kewley & Dopita 
l2002h . The variations of [N Il]/Ha a nd [N II]/[Q III] presen ted here 
show a similar trend as those given bv lErcolano et al J ( 120071) for dif- 
ferent configuration of stars and gas in H II re gions at the metallicit y 
of NGC 595 (Z= (0.4-1) Z Q i n the models oflErcolano et al.l2007b . 
Following the calibrations of IPettini & Pagell j2004h for [N II]/Hq 
and [Nll]/[OlIl] we find variations of ~0.3 in log[0/H] for both 



emission-line ratios. Finally, [Nll]/[Sll] and [Nll]/[Oll] are more 
stable within the region; in particular [Nll]/[Oll] would show a 
rather similar behaviour as R23, which would make it as good 
emission-line ratio to trace the metallicity of the regions as R23 
is. 



5 CONCLUSIONS 

We present IFS of NGC 595, one of the most luminous Hll re- 
gions in M33, covering an unprecedented area in the disc of the 
spiral galaxies of the Local Group, a ~ 174x340 parsec 2 field of 
view representing the complete surface of the region. Taking ad- 
vantage of the power of these observations, we are able to identify 
and catalogue the WR population of the region and to make the best 
census of WR stars in NGC 595 up till now. We have analysed the 
variations within the region of the main emission-line ratios that 
describe the physical properties of the region. The analysis of the 
observations yields the following results. 

• We present the properties of NGC 595 derived from the inte- 
grated spectrum of the region. The electron density and the metal- 
licity are consistent with previously reported values in the literature 
using long-slit spectroscopy, covering just the most intense knots in 
the region. 

• The extinction map obtained from the Ha/H/3 emission-line 
ratio at each spaxel in the field of view presents a concentrated dis- 
tribution with the maximum located at the centre of the Ha shell 
structure. We have compared this map with the Spitzer 24 and 8 /im 
bands using elliptical radial profiles. The 24 /im emission and the 
Ha/H/3 radial profiles are similar, with their maxima located at the 
same distance from the ionizing stars. The maximum of the 8 fim 
emission radial profile shows a displacement of ~4-5" (15-18 pc) 
with respect to the maximum of the Ha/H/3 radial profile. This 
shows that the extinction suffered by the gas is produced by dust 
emitting at 24 /im, which is probably mixed with the ionized gas 
within the region, while the dust emitting at 8 ^im is more related to 
the outer parts of the region where the molecular cloud is located. 

• The [S II]A6717/[S II] A673 1 emission-line ratio map does not 
show any structure, implying that the electron density is quite con- 
stant within the region, despite the pronounced Ha shell morphol- 
ogy of NGC 595. The values of [S n]A6717/[S n]A6731 are consis- 
tent with the low electron density limit. 

• We have produced BPT diagrams using the emission-line ra- 
tios at each spaxel in the field of view. We find that the spaxels 
with a low Ha flux cover the total excitation range of the region. 
The location of the spaxels with a high Ha flux in the [O IIlJ/H/3- 
[S Il]/Ha diagram is displaced with respect to the position given 
in the diagram for the integrated values. This shows the limitations 
of the long-slit observations, which usually cover the most intense 
Ha knot in the regions, to obtain a reliable value for the [S Il]/Ha 
emission-line ratio. 

• The study of the WR population for NGC 595 has been com- 
pleted with these observations. We have covered the total surface 
of the region and found a new WR candidate whose spectrum is 
characteristic of a WN star. The WR candidate, situated far away 
from the other WR stars in the region, is close to the Ha shell but 
corresponds to a zone of low Ha emission. 

• We have analysed the behavior of the most used metallicity 
calibrators: R23 shows small variations of 0.13 dex despite the spa- 
tial distribution of the stars and gas in the region and the strong 
trend of [O IIl]/[0 II] to decrease to the outer parts of the region. 
This shows the robustness of R23 to estimate the metallicity of 
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Figure 12. Top: maps of the oxygen emission-line ratios ([O III] AA4959, 50()7/[O II] A3727) tracing the ionization parameter (left-hand panel) and the 
metallicity ( R23= ([O II] A3727 + [O III] AA4959,5007)/H/3) (right-hand panel). The emission-line ratios shown here are those having relative errors < 30 per 
cent. The orientation is the same as in Fig. [6] Bottom: elliptical radial profiles of [OlII]/[OlI] (left-hand panel) and R23 (right-hand panel, black diamonds). 
The ellipse parameters are the same as the one used in Figs.[5]and |7]and the rings have 2 arcsec width. In the right-hand panel, we also plot radial profiles of 
the main emission-line ratios used as metallicity calibrators: [N II]/[0 II] (blue), [N II]/[S II] (green), [N II]/Hcj (red) and [N II]/[0 III] (magenta). 
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the region; in spite of the H II region morphology, far away from 
the classical Stromgren sphere picture, the R23 parameter varies 
slightly within the region and thus can be used to obtain a rep- 
resentative value of the H II region metallicity. In contrast, other 
parameters such as [N Il]/Ha and [N Il]/[0 III] show strong varia- 
tions within the region (up to one order of magnitude in the case of 
[Nll]/[0lll]), which show their strong dependence on the ioniza- 
tion parameter. 

We show in this paper the power of the IFS, which overcomes 
the limitations of long-slit observations in H II regions. The results 
presented here show that for the emission lines dominated by the 
contribution of the most intense knots, long-slit observations are 
representative of the integrated emission of H II regions. The situa- 
tion is different in the case of lower excitation emission lines, as it 
happens in the [S II] emission lines, for which there is a significant 
bias for long-slit observations. This bias has only been possible to 
quantify using IFS observations with a complete H II region cover- 
age. 

The physical properties can vary within H II regions and the 
emission-line ratio maps derived from IFS observations allow us 
to make a detailed study of these variations. For NGC 595 we find 
that the ionization parameter shows strong variations with the radial 
distance from the stars; the ionization structure is clearly depicted 
with the corresponding emission-line diagnostics and the reddening 
map presents a non-uniform distribution with a maximum located 
at the center of the observed Ha shell structure. The electron den- 
sity map, however, does not present significant variations within the 
region. 

A change of one order of magnitude in [Oll]/[OlIl] within 
NGC 595 allows us to study the dependence of the ionization pa- 
rameter on the most widely used metallicity calibrators, since other 
properties such as density or the ionizing spectrum do not change 
within this H II region. The main result is that R23 does not depend 
significantly on the ionization parameter, but other emission-line 
ratios such as [N II]/Hq and [N Il]/[0 III] show important variations 
within the region. 

Finally, the capability of the IFS data to make complete cen- 
suses of the WR population in stellar clusters is very nicely shown 
here: while the hunting of WR stars are based on identification in 
broad-band images with follow-up spectroscopic observations, we 
are able to perform such an analysis in only one night of obser- 
vation, covering the whole surface of the region and allowing us 
to identify WR stars that are not located close to the centre of the 
cluster. 
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